library(knitr)
library(tidyverse)
library(lifecontingencies)
library(actuar)

#导入《中国人身保险业经验生命表（2010－2013）》
#资料来源：http://www.isc-org.cn/jgzc/5237.jhtml
clt<-read.csv("lifetable/lifetable_cn.csv",
              header=T)
knitr::kable(head(clt,6))

ggplot(data=clt)+
  geom_line(aes(x=age,y=CL1,color="非养老1：男"))+
  geom_line(aes(x=age,y=CL2,color="非养老1：女"))+
  geom_line(aes(x=age,y=CL3,color="非养老2：男"))+
  geom_line(aes(x=age,y=CL4,color="非养老2：女"))+
  geom_line(aes(x=age,y=CL5,color="养老：男"))+
  geom_line(aes(x=age,y=CL6,color="养老：女"))+
  xlab("年龄")+
  ylab("死亡率")+
  ggtitle("中国人身保险业经验生命表（2010－2013）")+
  theme(legend.position = c(0.15,0.55)) 
ggplot(data=clt)+
  geom_line(aes(x=age,y=CL1,color="非养老1：男"))+
  geom_line(aes(x=age,y=CL2,color="非养老1：女"))+
  geom_line(aes(x=age,y=CL3,color="非养老2：男"))+
  geom_line(aes(x=age,y=CL4,color="非养老2：女"))+
  geom_line(aes(x=age,y=CL5,color="养老：男"))+
  geom_line(aes(x=age,y=CL6,color="养老：女"))+
  xlab("年龄")+
  ylab("死亡率")+
  ggtitle("中国人身保险业经验生命表（2010－2013）")+
  theme(legend.position = c(0.15,0.55)) 
# ggsave("lifetable/lifetable.png",width=400,hight=300,units = "px")


cl1 <- clt %>% select(age,CL1)
cl1 <- probs2lifetable(clt$CL1,radix=10^5,type="qx",name="lt_CN")

print(cl1)
head(cl1,3)
cl1@lx[1:4]
cl1@x[1:4]

?exn
